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This is a study of a monochromatic planar perturbation impinging upon a canonical acoustic 
hole. We show that acoustic hole scattering shares key features with black hole scattering. The 
interference of wavefronts passing in opposite senses around the hole creates regular oscillations 
in the scattered intensity. We examine this effect by applying a partial wave method to compute 
the differential scattering cross section for a range of incident wavelengths. We demonstrate the 
existence of a scattering peak in the backward direction, known as the glory. We show that the 
glory created by the canonical acoustic hole is approximately 170 times less intense than the glory 
created by the Schwarzschild black hole, for equivalent horizon-to-wavelength ratios. We hope that 
direct experimental observations of such effects may be possible in the near future. 

PACS numbers: 04.70.-s, 11.80.-m, 47.35.Rs 



I. INTRODUCTION 

Black holes play an important role in modern physics. 
For many years, they were considered merely a mathe- 
matical curiosity arising from Einstein's field equations. 
In more recent decades there has emerged strong evi- 
dence for their presence in binary systems and galactic 
centers. Aside from the growing observational data, it is 
clear that the idea of a black hole has its own attractive 
power. Black holes have captured the imagination of a 
generation of physicists. For example, debates on black 
hole thermodynamics [T] and the information paradox [5] 
are ongoing. The detection of gravitational waves from 
black hole binaries is hoped to be imminent [3,. The 
possibility of black hole creation at the LHC [H [5] has 
recently received hyperbolic media attention. 

Naturally, the possibility that black hole analogues [,6j 
may be created in the laboratory has attracted 
widespread interest. Analogues are artificial systems 
with some of the key properties of black holes. In 1981, 
Unruh [7j proposed the idea of a dumb hole: a region of 
fiuid from which no sound may escape. A dumb hole, or 
acoustic hole [S], is bounded by an apparent horizon: a 
surface through which the normal flow velocity is equal 
to the speed of sound in the fluid. In this paper, we 
study the simplest spherically-symmetric dumb hole, the 
so-called canonical acoustic hole. A range of alternative 
analogue systems have been proposed, for example in su- 
perfluid helium [S] , in Bose-Einstein condensates [TD] , in 
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electromagnetic waveguides pjj , in optical fibers [T^l [TB] 
and other systems [T3]. 

Under certain assumptions (an inviscid, barotropic 
fiuid; irrotational fiow), the Navier-Stokes equations de- 
scribing small perturbations to fiuid fiow turn out to be 
formally identical to the equations for a massless scalar 
field propagating in a (3-1-1) Lorentzian geometry. That 
is [8 , 

□ ^ = -^a^ {^ggt^-d,^) = 0, (1) 

where ip is a. small perturbation to the potential de- 
scribing the fiuid velocity, v = Vq — V^". The effective 
metric g^j/ is an algebraic function of the background 
fiow Vq, the local density p{P), and the speed of sound 
— dP/dp, where P is the fiuid pressure. The appro- 
priate metric for the canonical acoustic hole is given in 
section |TT1 

Acoustic black holes are expected to emit a ther- 
mal spectrum of phonons, since Hawking's original argu- 
ments follow through without essential modification [7]. 
On the other hand, the laws of black hole thermo- 
dynamics do not follow. It seems that the associa- 
tion of entropy with horizon area is a special prop- 
erty of the Einstein equations [T^. The Hawking 
temperature for the canonical acoustic hole, Tu « 
1.2 X lO-^K m(c/1000ms-i)(c-idu/dn) |8j, is suffi- 
ciently small that direct detection of Hawking radia- 
tion has not been achieved. The possibility of detecting 
Hawking emission from a ID Bose-Einstein analogue is 
under investigation [TB] . 

There is more to acoustic holes than just Hawking ra- 
diation, however. As has been pointed out |17j . there 
is a range of interesting classical wave phenomena which 
should be amenable to experiment. In particular, the 
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scattering of perturbations by acoustic holes has been 
studied by a number of authors, in both the frequency 
domain [THl UHl 120] and the time domain [5T]. Like 
black holes, acoustic holes have characteristic damped 
resonances called quasi-normal modes. These resonances 
appear in the dynamic response of acoustic holes to ex- 
ternal perturbation. Time-dependent scattering has been 
considered in detail for the lowest multipoles / = 0, 1, . . .. 
However, a study of the overall scattering pattern pro- 
duced by the superposition of many angular modes has 
not been considered. This is the main purpose of our 
paper. 

In this study, we adopt the partial wave approach to 
investigate the scattering and absorption of a monochro- 
matic external perturbation impinging upon a canoni- 
cal acoustic hole. We compare our results with previous 
studies of time-independent scattering by astrophysical 
holes [25]. In this regard, we are motivated by the excit- 
ing prospects for gravitational wave astronomy [3]. It is 
possible that acoustic scattering experiments may offer a 
way to better understand the dynamics of gravitational 
radiation. Further studies of wave propagation in ana- 
logue systems and their experimental realization are to 
be encouraged. 

The monochromatic scattering scenario depends on 
only a single dimensionless parameter |23j , 



(2) 



which expresses the ratio of the incident wavelength A to 
the horizon size rh of the canonical acoustic hole. For 
convenience, we have set the speed of sound c equal to 
unity from here on. 

The remainder of the paper is organized as follows. 
In section |TT] we introduce the canonical acoustic hole, 
and the equations which describe the propagation of 
sound perturbations in the effective geometry. In sec- 
tion |TTT] we examine the properties of geodesies on a class 
of spherically-symmetric geometries. We derive strong 
and weak-field approximations for the scattering angle, 
and extend the glory approximation [24 to the case of 
the canonical acoustic hole. In section IIVI we review 
the partial-wave method for time-independent scattering. 
An approximation for the phase shift in the large-Z limit 
is obtained. In section |V] we discuss the details of the 
numerical methods employed. A selection of key results 
is presented in sec tion [yi] We conclude with our final 

remarks in section VII The metric signature (H ) 

is used throughout. 



II. THE CANONICAL ACOUSTIC HOLE 

Let us consider a spherically-symmetric steady flow of 
incompressible fluid in three dimensions, with a source 
or sink at the origin r = 0. Conservation of fluid implies 
a radial velocity = ±r^/r^, where r/j is the radius at 
which the flow speed exceeds the speed of sound in the 



fluid. Perturbations 5v = — ViA to the steady flow are 
governed by Eq. ([T]) with the effective geometry 



ds^ = g^.^dxt'dx" = +dP - [dr"^ ± ^dtj 

-r^ {de^ + sin^ 0d(P'^) . (3) 



The {+) sign stands for a source and the (— ) sign stands 
for a sink. We will henceforth choose the (— ) sign, so that 
the system is analogous to a black (rather than white) 
hole. By introducing an alternative time coordinate dt = 
dt — {rf^/r'^){i — r'^/r'^)~^dr the metric may be written 
in a diagonal form. 



ds^ 



+f{r)df - r\r)dr^ - r^{de^ + sin^ 0d(j)^), (4) 



where f{r) = 1 — r^/r'^- Note that the time parameter t 
diverges as r ^ rh- 

Since the perturbation equation ([ij is simply the 
Klein-Gordon equation on a curved background, we may 
employ standard field-theory techniques. We allow the 
field -0 to take complex values, and construct the con- 
served current 



satisfying 



{g^'^J,) = 



(5) 



(6) 



Alternatively, for direct comparison with acoustic exper- 
iments we may work directly with the real part of the 
field Re(?/j), which corresponds to perturbations in the 
flow velocity, (5v = — V[Re(V')]. 

The appropriate choice of observable ( or -0) depends 
on the point of view adopted. If we consider acoustic 
scattering experiments to be a model for black holes irra- 
diated by incoherent electromagnetic radiation, it makes 
sense to work with J^, which behaves like an intensity. 
If, on the other hand, we are interested in coherent grav- 
itational radiation impinging on a black hole, we should 
consider the fleld itself. Interferometers such as LIGO 
will measure a dimensionless strain h = l^L/L ~ l/r 
rather than an intensity / ~ 1 /r^ |25i . In the remainder 
of the paper, we will take the former approach, and apply 
standard techniques from quantum mechanics |26j . 



III. GEODESICS AND ORBITS 

In the short-wavelength limit (A rh), perturbations 
obey the eikonal approximation ip ^ e»'=''^f< ^ where fc'' is 
tangent to a null geodesic {k^k^ = 0) on the background 
geometry. This naturally motivates a study of the scat- 
tering of null geodesies on spherically-symmetric, static 
backgrounds. For greater insight, let us keep the analysis 
general and consider the class of geometries described by 
line elements 



ds' 



-fdt^ - r^dr^ - r^{de^ + sin^ 



(7) 
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where / = 1 — r^/r"^ and Vh is the horizon radius of the 
hole. Here, n = \ corresponds to the Schwarzschild black 
hole (a solution of the vacuum Einstein field equations) 
and n = 4 corresponds to the canonical acoustic hole (a 
solution of fluid-flow equations) . 

The paths of scattered geodesies are found from the 
solutions of the orbital equation 



du 



1 
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(8) 



where u = 1/r, and b is the impact parameter. Differen- 
tiating Eq. ^ leads to a version of Binet's equation [27] . 



n n+l 



(9) 



Solutions of Eq. (|8| give the deflection (or scattering) 
angle as a function of the impact parameter, 9{h). For 
certain values of n, closed- form solutions may be found in 
terms of elliptic integrals. For example, for the canonical 
acoustic hole (ti = 4), the scattering angle is 



e{b) 



2K{k) 

y/{v2 - Vl)v3 



- 71", 



(10) 



where 



W2(U3 - vi) 



V3{V2 - Vi) ' 

Here, Vi, V2 and are the roots of the cubic 
rf,v^ -v + l/b^ = 0, 



(11) 



and K{k) is a complete elliptic integral of the first 
kind [28 . For scattering trajectories, we have vi < 
0, > and > V2- 

Let us consider approximate solutions of Eq. ([8| in the 
weak- and strong-field regimes. By treating the term on 
the right-hand side of Eq. (|9| as a small perturbation, it 
is straightforward to show that, in the weak-field limit, 
the defiection angle 9 is approximately 



(12) 



where the constant of proportionality is a„ = 
{2,37r/4,8/3,157r/16,16/5,357r/32} for n = 1...6, re- 
spectively. In particular, for the Schwarzschild black hole 
(n — 1) we get Einstein's defiection angle 9 « ^r^/b, 
and for the canonical acoustic hole {n = 4) we find 

An unstable circular orbit is present at the critical ra- 
dius Tc — (^^^y^) f*/!. This in turn implies a critical 

^ /2 

impact parameter, be — Vc {I — r^i^/r'c) ■ An incident 
ray starting with impact parameter b = be ends on the 
unstable orbit &t r = r^, the ray circles the hole an in- 
finite number of times. Rays with b < be are absorbed 
by the hole; rays with b > be are scattered. The classical 



n 


Vc/rh 


bc/rh 




1 


1.5000 


2.5981 


6.7500 


2 


1.4142 


2.0000 


4.0000 


3 


1.3572 


1.7522 


3.0700 


4 


1.3161 


1.6119 


2.5981 


5 


1.2847 


1.5201 


2.3108 


6 


1.2599 


1.4548 


2.1165 



TABLE I: Numerical values of the unstable orbit radius (vc), 
the critical impact parameter (6c) and the absorption cross 
section (aabs) for n = 1 ... 6, where n — 1 corresponds to 
the Schwarzschild black hole, and n = 4 corresponds to the 
canonical acoustic hole. 



absorption cross section is simply the area of a circle of 
radius be, that is dabs — ^fe^. Numerical values for be 
and (Tabs are listed in Table [ij 

Geodesies passing very close to the unstable circular 
orbit may be scattered through large angles. For b> be, 
the defiection angle is approximately 



(13) 



Here Pn is a numerical coefficient which may be calcu- 
lated on a case- by-case basis. Many years ago, Dar- 
win [29] showed that /3i « 0.6702 for the Schwarzschild 
black hole. By considering the asymptotic form of 
Eq. (lOl as fc ^ 1, it is straightforward to show that 

= 108e^2-n' ~ 0.2017 for n = 4 (the canonical acoustic 
hole). The calculation is outlined in Appendix [A| 

It is well-known that an unstable orbit gives rise to a 
glory in the backward direction. A glory is a bright spot 
or ring in the backward scattering direction. The size and 
intensity of the spot or ring depends on the wavelength. 
Via path-integral arguments it was shown that ^3] , close 
to 9 — TT, the scattering cross section is approximately 



da 

dn 



2nujbl 



[J2s{bgUJ sin ( 



(14) 



Here bg is the impact parameter for which the defiection 
angle is tt, J2s{x) is a Bessel function, and s is the spin 
of the scattered field (in this case s = 0). 

Combining results ([T3| and ( 14 1 leads to the glory ap- 
proximation 



, da 

dn 



0.019835 Lur^ [ Jo(1.61246 ujr,, sin 9)f (15) 



for the canonical acoustic hole [30 , which may be com- 
pared with the approximation 



rii 



, da 

dh 



3.3772 Lurh [ Jo(2.67325 ur,, sin 9)f (16) 



for the Schwarzschild black hole. In other words, the 
glory peak for the canonical acoustic hole is similar in 
nature to the Schwarzschild black hole, but is approxi- 
mately 170 times weaker in magnitude |31j . For a given 
frequency, the angular width of the glory peak of the 
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Schwarzschild black hole is narrower than for the canon- 
ical acoustic hole, by a factor of approximately 1.66. 
In the next section, we test the validity of approxima- 
tion ( 15 1 by computing the exact cross section via partial 



wave series. 



IV. PARTIAL WAVE METHOD 

In this section we use the partial wave method to de- 
termine the differential scattering cross section da/<Kl at 
intermediate wavelengths (A ^ r^). 

By substituting the metric Q and the separation 
ansatz 



r 

into Eq. ([l| we obtain a radial equation 
where 

7' , 1(1+ 



Vi{r) = f[- + 



(17) 



(18) 



(19) 



and /' = df /dr. The perturbation must be purely ingo- 
ing at the horizon r — r^- Hence we impose the boundary 
condition 



u{r) = exjp(—iujx) 



(20) 



where 



so that 



dx 
dr 



rii 



In 



Th tan 

2 



4 ' 



is the Wheeler-type coordinate. The constant of integra- 
tion has been chosen so that x « r as r — > cx3. It is 
straightforward to check that in the original coordinate 
system ^ the ingoing solution e~*'^^ is well-defined at 
r = Th, whereas the outgoing solution e''""^^ is divergent. 
Towards spatial infinity, the asymptotic form of the so- 
lution is 

ui{r) ~ ujx [Aout i^+^h\'\uox) + A,^ {-iY^+^^h^^^* {lox) 

(21) 

where hf^{x) are the spherical Bessel functions of the 
third kind |28j . and A^^ and Aout are complex constants. 
The phase shifts 5i are defined by 



(-l)'+MoutMi, 



(22) 



and the scattering amplitude is 

2i^E(2/-Hl)(e2^^'-l)P,(cos( 



(23) 



From the amplitude, the differential scattering cross sec- 
tion follows immediately. 



(24) 



Finally, we can define the scattering (elastic), absorp- 
tion (inelastic) and total (combined) cross sections [25] . 



'^scat 



da 



dn 

oo 
1=0 



-dn 



(21+1) \e 



2i5i 



fabs 



CTtot 



oo 

JB2'+1) (l-|e^"f) 



1=0 

f^scat + *^abs 

27r 



(25) 
(26) 



UJ 



J2i2l+l){l-Re{e'^'^)) . (27) 



To compute the cross sections (24 27) we must 



solve Eq. lYiM subject to the boundary conditions (20 1 



and (21 1, to obtain numerical values for the phase shifts 



via Eq. ( 22 1 . The numerical method employed is de- 



scribed in section |V] The task is made substantially eas- 
ier if we utilise an approximation for the phase in the 
large-/ limit. We find that, for Z /c ^ ^^bc, the phase 
shift for the scattering from the canonical acoustic hole 
is approximately 



Si 



32(/+ 1/2)3 ■ 



(28) 



There are (at least) two ways to reach this result. 

The first way to obtain Eq. ( 28 1 is through Ford 
and Wheeler's [32] semiclassical description of scatter- 
ing, whereby an impact parameter is associated with each 
partial wave. 



6 = 



/+ 1/2 



(29) 



Much physical information can be extracted from the de- 
flection function defined by 



6(0 = - [Re{2Si) 



(30) 



Here, / is allowed to assume continuous real values. 



Equating the deflection function ( 30 1 with the weak- 
field deflection angle {12) yields 



dl 



[M26i 



157r 157r (wr/j)^ 

l6~F " "T6~(/ + 1/2)^ 



(31) 
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from which Eq. ( 28 ) follows immediately (up to a sign) . solution onto Eq. ( 20 1 . We checked the accuracy of the 



The second way to obtain Eq. (28) is through the 
Born approximation First, we make the substitu- 



tion u{r) = f~^/'^X{r) in Eq. so that 



+ r 



fd'fl df 



2 dr^ 4 V dr 



X ^0. 



(32) 

Next, we approximate the wave equation as a power series 
in 1/r, 



d^X 
dr^ 



LO^ - ^ „ ' + U{r) 



X = 0, 



(33) 



where 



o 



,10 



Now we apply the Born approximation formula 



5i « —Lo I r ii {ijjr)U{r)dr, 
Ja 



(34) 



(35) 



where ji{x) are the spherical Bessel functions of the first 
kind [2H] . The first two terms in U (r) give a contribution 
proportional to {ujrh)'^. Applying the identities 



(2?-l)(2;+l)(2/ + 3) 

oo 



X ii {x)dx ■ 







371- 



{21 - 3)(2/ - 1)(2Z + 1)(2/ + i){2l + 5) 



(36) 



(37) 



leads to Eq. (28 1 to lowest order in / + 1/2. 



NUMERICAL METHOD 



To compute numerical phase shifts, we first solve the 



radial equation (18) subject to Eq. (20) and compute 



ingoing and outgoing coefficients {A-^-^ and Aouti respec- 
tively) by matching on to Eq. ( pTj ). The phase shifts 
follow immediately from Eq. ( 22 1 . 



Our numerical method is straightforward. We start 
close to the horizon, at r — + rjrh (where typically 
r] ~ 0.001) with a series expansion 



k=0 



(38) 



The coefficients ak can be found with a symbolic algebra 
package such as Maple or Mathematica. We evaluated 
the series at 6th order {km = 6). Next, we numerically 
integrate from near the horizon into the asymptotically 
flat region, r = '--^ ISOr^j. We matched the numerical 



resulting phase shifts by varying rj and r™. 

Fig. [T] shows typical phase shifts as a function of angu- 
lar mode number I. The real and imaginary components 



of 



are plotted. It is clear that partial waves of low 



order {I ^ 1^ — Lobc) are almost entirely absorbed. Par- 
tial waves of intermediate order {I ~ Ic) are partially 
scattered and absorbed. In this regime the phase shift 
varies rapidly with I. For partial waves of large order 
{I ^ Ic) the phase shift 6i is purely real, and conforms to 
approximation ( 28 1 . 



Phase shifts for a canonical acoustic hoie at ti)r(|=6.0 

I I I I I I I I I 




mode number 1 



FIG. 1: Phase shifts for a canonical acoustic hole at Lor^ = 
6.0. Below Ic ~ Ldbc ~ l.Q12ujrh ~ 9.7, absorption dominates. 



In the large-/ limit, e 



l + ^iTi{urhY/l&{l + l/2Y 



In practice, we must truncate the partial wave sum 
in Eq. ( 23 1 at some finite upper limit / = /,„ . If done 



naively, this introduces a truncation error. To avoid this 
problem, we employed the asymptotic phase shifts (28) 



to minimize the loss of precision. The amplitude may be 
split into two parts, 

f{0) = .ftrli0)+fj.i^{9) 



— Y{2l+l){e^'^'~l)Pi{cos( 

if, J ^ — ^ 

5Tru;^rf;\ ^ Pi{cos0) 



2iLU 



16 



^ (^+ 1/2)2 



(39) 



The remainder term frem may either be approximated 
numerically (by choosing some extremely large cutoff) 
or, in certain cases, computed analytically. For example, 
to compute the sum in the forward direction 9^0, for 
which P;(l) = 1, we may use the result 



J2 {1 + 1/2)-^ = ^{l,lm + 3/2), 



(40) 



1=1^+1 



where ^! is the polygamma function f28J. To compute 
the remainder sum in the backward direction 9 — tt, for 



6 



which Pi{~l) = (—1)' we may use the result 

oo 

^(-l)'G + l/2)-2 = 4C, (41) 



Scattering cross sections at Djr^=6.0 



1=0 



where C ~ 0.91596559 is Catalan's constant [34]. 



VI. RESULTS 

In this section, we present our numerical results and 
compare with theoretical approximations and expecta- 
tions. 

Fig. |2]shows the scattering cross section as a function of 
angle, for a range of coupHngs (for ujrh = 0.2 to 6.0). At 
small couplings (wr/j <C 1), the scattered flux is isotropic, 
and the I = mode dominates. At higher couplings, 
flux is preferentially scattered in the forward direction, 
and a more complicated pattern arises. In particular, for 



> 1 we see an oscillatory pattern with an angular 
width inversely proportional to wr/j. Similar patterns 
are observed for black hole scattering ^21 [SB ISSJ 157] . 
An explanation for the physical origin of the oscillations 
can be found in [24j [38] . The oscillations arise from the 
interference of rays that pass in opposite senses around 
the hole (i.e. interference between rays scattered through 
9,2TT-e,27r + e, ... etc. [35]). 

Scattering cross sections for uir|, = 0.2 to 6.0 



scattering angle ft (deg) 



FIG. 2: Scattering pattern of a canonical acoustic hole for a 
range of incident wavelengths. The plot shows the (logarithm 
of the) differential scattering cross section da/dO, as a function 
of scattering angle 9 for a range of frequencies ujrh . 

Fig. |3] shows an example of the interference effects at 
ujTh = 6.0. The plot compares the scattering cross sec- 
tion for the Schwarzschild black hole with the canoni- 
cal acoustic hole. The Schwarzschild black hole scatters 
much more flux, due to the long-range nature of the in- 
teraction {(j) ~ 1/r versus ( 



1/r ). In addition, the 



oscillations from the Schwarzschild black hole are nar- 
rower than from the canonical acoustic hole, by a factor 
of approximately 1.66. The glory approximations (15) 



and (16 1 are shown as broken lines. 




Glory approximation 
Sctiwarzsctiild black tiole 
Canonical acoustic hole 



scattering angle 9 (deg) 



FIG. 3: Scattering cross section of Schwarzschild black hole 
and canonical acoustic hole at ujrh = 6.0. The numerically- 
determined scattering cross sections are compared with the 
glory approximations (Eqs. (15 1 and (16l). Note the loga- 



rithmic scale on the vertical axis. 



Fig. |4] shows the differential scattering cross section 
in the backward direction, 9 = 180°, as a function of 
frequency. At low frequencies, the back-scattered flux 
is predominantly due to the I — mode. At high fre- 
quencies, the glory approximation predicts that the flux 
increases linearly with wr^. The glory approximation is 
shown as a straight line, and provides a good fit to the 
numerical data. 



Flux scattered in ttie backward direction (fl = n) 




FIG. 4: Flux scattered in the backward direction. The plot 
shows da/dQ a.t 8 — 180° as a function of the coupling corh. 



We note that for the canonical acoustic hole, the scat- 
tering cross section in the forward direction {9 — 0°) is 
well-defined and finite (see Fig. |2]). This contrasts with 
the Schwarzschild black hole, for which da/dO. ~ 9~'^ as 
9 ^ Q. It also contrasts with the asymptotic behavior 
of the 'classical' scattering cross section, which may be 
defined with reference to parallel geodesies approaching 
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from infinity, 



Canonical acoustic hoie scattering cross section 



da 




b 


db 




cl. 


sin 9 


d0 



(42) 



Substituting in the weak-field deflection angle ( 12 1 leads 
to da/dn\,i, ~ 6'-^/2. 



Flux scatterect in the fonAfarcJ cjirection (H = 0) 




FIG. 5: Flux scattered in the forward direction. The plots 
show da/dfl at 6^ = 0° as a function of the coupling ujrh- The 
lower plot shows the same data as the upper plot, but with a 
log-log scale. The data suggests that the on-axis cross section 
diverges as (ujrh)'' with an exponent 3 < 7 < 4, as curh 00. 

The scattering cross section in the forward direction 
(9 = 0°) as a function of frequency is shown in Fig. [s] It 
increases rapidly and monotonically with wr/j. However, 
we should not forget that Eq. ([ij is only valid for small 
perturbations to the background flow. Hence it is likely 
the linear approximation will break down on axis. We 
would expect non-linear (second-order) effects to become 
important close to 9 = 0°. 

The overall scattering and absorption cross sections 
(defined in Eqs. ( psj ) and (26 1 ) are shown in Fig. |6]and 
in Fig. [7] respectively, for < ujrh < 6. As has been ob- 
served previously [23], the absorption cross section tends 
to the horizon area tXabs ~ 47rr^ at low frequencies |40j . 
and to the geometric-optics value a^hs ~ T^b^ at high fre- 
quencies. On the other hand, the scattering cross section 
increases with Lurh, in an approximately linear fashion. 



VII. FINAL REMARKS 

In this paper we have presented an investigation of the 
scattering of sound waves by canonical acoustic holes. We 
have applied analytic and numerical methods to deter- 
mine the differential scattering cross section for a range 
of wavelength-to-horizon ratios of a monochromatic pla- 
nar perturbation impinging on the canonical acoustic 



FIG. 6: Scattering cross section for < uirh < 6. The 
scatteri ng c ross section for the canonical acoustic hole, given 
by Eq. ( 25 1, is plotted as a function of the coupling corh- The 
contributions from the I = . . .7 partial waves are shown. 



Canonical acoustic hole 



ion cross section 



FIG. 7: Absorption cross section for < urn < 6. The 
canonical acoustic hole absorption cross section given by 
Eq. (261 is plotted as a function of the coupling urh- The 
contributions from / = up to 7 are also shown. 



hole. We obtained high-precision results for the scat- 
tering cross section with the partial wave method, con- 
sidering all angular momentum contributions. 

The diffraction pattern created by a canonical acoustic 
hole is distinctive (Fig. |2| . It shares some of the features 
of Schwarzschild black hole diffraction ^j. There are 
oscillations in the scattering amplitude (Fig. [3| due to 
the interference of wavefronts passing in opposite senses 
around the hole (so-called orbiting [38. ). The structure 
of these oscillations is intimately linked to the properties 
of the unstable null orbit which lies close to the horizon. 

Experimentally, such interference effects would be 
manifest as angular segments of alternately weak and 
strong intensity, whose angular width is proportional to 
wavelength. We expect a particularly strong segment 
of positive interference in the backward direction: the 
so-called glory. In this paper we computed the glory via 
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semi-classical approximations (Eq. ( 15 1) and partial wave 
methods (Fig. [3|, and found the two approaches to be 
consistent (FigTMl. 

Wave scattering by canonical acoustic holes is different 
to wave scattering by astrophysical black holes in three 
key respects. Firstly, at long wavelengths (ujrii ^ 1), 
scattering by the canonical acoustic hole is isotropic 
(Fig. |2]) whereas scattering by the Schwarzschild black 
hole is Coulombian, i.e. da/d^ « r^/4sin''(6'/2) [32]. 
Secondly, for equivalent couplings urh [31] , the intensity 
of the flux scattered by a canonical acoustic hole is much 
weaker than the intensity scattered by a Schwarzschild 
black hole, by approximately two orders of magnitude 
at large angles (Fig. |3|. Thirdly, the scattering ampli- 
tude for the canonical acoustic hole is finite in the for- 
ward direction (Fig. |2]), whereas the amplitude for the 
Schwarzschild black hole diverges as 6* — > [22 . These 
differences are explained by the behavior of the effec- 
tive potential ( 19 ) far from the hole. The potential of 
the canonical acoustic hole decays more rapidly (^ 1/f^) 
than the 'Newtonian' potential (~ I/?')- 

Given the exciting prospects for gravitational wave as- 
tronomy [3] , time-domain simulations of wave scattering 
by acoustic holes would be of interest . Even better would 
be an experimental realization of wave scattering by the 
canonical acoustic hole, or indeed, by any of the alterna- 
tive black hole analogues [SlITKHliaiiniinilllllKIll- 
We propose this as a challenge for experimental groups 
worldwide. 
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APPENDIX A: LOGARITHMIC SCATTERING 
APPROXIMATION 

Here the constant /3„ in the spiral scattering approxi- 
mation (13 1 is computed for the special case n = 4 (the 
canonical acoustic hole). We start with the exact so- 
luti on ( 10 1. The roots wi,W2,f3 are defined by the cu- 
bic (111. The critical orbit occurs when rf^vi = — 2/\/3, 
and rlv2 = rlva = l/VS, and b = = {3^/'^ /2^/^)rh. 



For orbits close to the critical orbit, the roots are 



rlvi 
rlv2 
rlv3 



I/31/2 - 25/43-"/«<5 + 0{S^) (Al) 

1/31/2 + 25/43-11/8,5 + 0{d^). 



where S = [{b — bc)/rh]^^^ . Applying the approximation 



K{k) w In 



(l-fc2)l/2 



(A2) 
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to Eq. (10 1, it is straightforward to show that 



1 



■In 



where /34 = lOSe-^'^ « 0.201684 



b-br 



f3ibc 



(A3) 
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